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This  report  preseats  eqnstioas  which  describe  the  deforaatioa  sad  aeckiag 
io  s  rod  stretohiag  uader  s  large  velocity  gradieat.  The  equations  duplicate 
wave  propagation  computer  program  solutions  for  rods  of  varying  densities* 
yield  strengths,  and  other  properties  being  stretched  under  uniform  velocity 
gradients. 

The  equations  provide  insight  into  similar  phenomena  seen  in  stretching 
jets  from  shaped  charge  munitions  and  provide  a  purely  mechanistic  approach  to 
understanding  such  phenomena.  They  suggest  tailoring  of  jet  gradients  and 
material  properties  to  achieve  specific  weapon  objectives. 

The  basic  equations  are  solved  for  right  circular  cylinders  stretching 
under  a  linear  velocity  gradient.  Nonlinear  gradients  and  shapes  other  than 
constant-diameter  rods  are  modelled  by  the  basic  differential  equations. 
General  solutions  are  not  presented  for  these  eases  since  the  equations  must 
be  solved  numerically. 


ncnoM  II 

■OD  SISBIOUIB  CALCOLmCNS 


Hie  TOOOY  two-dlaeiiBloiuil  Lagrengien  Wave  Propagation  coapnter  prograa 
(Reference  1)  was  need  to  calculate  deforaation  In  stretching  rods  of  initial¬ 
ly  constant  dlaaeters.  Several  calculations  were  undertaken  with  varying 
gradient,  length,  density,  yield  strength,  and  sound  speeds.  Several  of  these 
calculations  vill  be  discussed  In  this  section  to  provide  an  introduction  to 
the  basic  phenoaenology. 

Figure  1  shows  the  initial  grid  structure  for  a  lO-ca-long,  l-ca-diaaeter 
rod.  The  top  figure  is  drawn  with  the  saae  scale  on  both  axes.  The  z,  or 
radial,  direction  is  aaplified  in  the  bottoa  figure.  The  z  axis  is  an  axis  of 
rotational  syaaetry. 

In  this  calculation  the  rod  is  copper  with  a  density  of  8.9  gn/cc,  an 
initial  bulk  sound  speed  of  4x10*  on/ sec,  a  shock  velocity/particle  velocity 
slope  of  l.S,  a  Qruneisen  ratio  of  2.0,  a  Poisson's  ratio  of  0.35,  and  a  Von 
Rises  yield  strength  of  5  kilobars  (Kb).  The  yield  strength  is  not  allowed  to 
vary  with  strain  or  internal  energy.  At  startup  tine  a  velocity  gradient  of 
2x10*  sec  ^  is  imposed  on  the  10-cm-long  rod.  The  velocity  is  set  to  0  at  the 
X>0  point  and  varies  linearly  with  X  to  2x10*  cm/ sec  at  the  X-10  cm  point. 

Figure  2  shows  the  rod  in  both  the  e^ual  and  amplified  scale  along  with 
X-velooity,  in  cm/ sec,  versus  X,  in  cm,  at  20  microseconds.  The  plot  clearly 
shows  relief  waves  moving  into  the  rod  from  both  ends,  stabiliiing  velocity  in 
the  relieved  sections.  The  velocity  in  the  relieved,  or  elastio,  sections 
increases  from  initial  values  in  the  tail  (low  velocity  end)  of  the  rod  and 
decreases  from  initial  values  in  the  tip  (high  velocity  end)  of  the  rod. 
Between  these  constant  velocity  segments  the  rod  is  continuing  to  stretch 
plastically  at  the  same  velocities  initially  given  to  the  rod. 
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Figure  3  presents  the  same  data  for  the  rod  at  60  microseconds.  At  this 
time  the  stabilized  end  sections  have  grown  to  approximately  2  cm,  with  a 
slight  neck  visible  at  the  end  of  the  relieved  sections.  The  velocity  and 
amplified  grid  plot  show  another  set  of  necks  farming  approximately  2  cm  ahead 
of  the  first  necks. 

By  100  microseconds  (Figure  4),  the  second  necks  are  well  formed  and  two 
new  necks  are  beginning  to  form.  The  stabilized  elastic  end  sections  are  each 
approximately  2.5  cm  in  length  and  an  elastic  region  also  exists  within  the 
adjacent  segment,  as  signified  by  the  constant  velocity  region.  Figure  5 
plots  axial  stress  in  dynes/cm*  versus  axial  position  along  the  rod  for  the 
zones  which  exist  in  the  radial  direction  across  the  rod. 

In  the  central  rod  section  stress  is  at  the  SxlO*  dynes/cm*  (5  Cb)  level. 
It  rises  above  this  level  within  each  neck  region  and  decays  below  the  level 
in  each  elastic  region.  The  rise  it  caused  by  an  increase  in  pressure  (mean 
stress)  due  to  the  curvature  at  the  neck.  This  rise  in  mean  stress  is  similar 
to  that  measured  in  tensile  tests  and  described  by  P.  V.  Bridgman  in  the  early 
1950's  (Reference  2).  As  discussed  in  Section  III,  it  is  this  rise  in  mean 
stress  which  causes  more  than  the  first  necks  to  form.  Figures  6  and  7  are 
plots  of  axial  strain  and  internal  energy  per  unit  mass  versus  axial  position 
at  the  100-microsecond  time.  Strain  along  the  completely  plastic  central 
section  of  the  rod  is  predictable  from  the  gradient  and  is  1.1  (110  percent). 
Axial  strain  rises  at  each  neck  location  and  decays  in  each  elastic  segment, 
as  would  be  expected.  The  strain  level  at  the  very  ends  of  the  rods  is 
virtually  zero  since  these  ends  were  relieved  almost  instantaneoualy.  The 
internal  energy  density  plot  in  Figure  7  shows  similar  phenomenology.  The 
energy  is  a  constant,  predictable  value  along  the  unrelieved  central  rod 
section  and  rises  as  strain  rises  in  each  neck.  In  this  calculation,  the  rod 
began  at  zero  internal  energy  at  the  start  of  the  problem. 

Figure  8  plots  three  times  the  square  root  of  the  second  stress  invariant 
(i.e.,  the  quantity  which  cannot  exceed  the  Von  Hlses  yield  strength)  in 
dynea/cm*  versus  axial  position  at  100  microseconds.  It  demonstrates  that  the 
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Figure  3.  Grid  and  Velocity  Plot  at  60  Microseconds 
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Figure  4.  Grid  and  Velocity  Plot  at  100  Microseconds 
Copper  with  2x10^  sec  ^  Gradient 
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Figure  6. 
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rod  is  Indsed  olsstie  in  the  end  sections,  plastic  at  the  necks  and  the  still¬ 
stretching  central  rod  section  and  elastic  in  the  second  set  of  stabilized 
segments  from  each  end.  Variations  at  the  sane  axial  position  are  Tariations 
within  zones  at  different  radii  (Z  position)  at  that  axial  position. 

Fignre  9  is  a  plot  of  kinetic  energy  density  (kinetic  energy  per  unit 
mass)  wersns  axial  position  at  100  microseconds.  It  shows  again  the 
stabilized  velocity  sections  and  the  cnrvatnre  in  the  center  section  expected 
from  a  constant  velocity  gradient. 

Necks  continne  to  multiply,  moving  in  from  each  end,  until  all  of  the  rod 
has  been  subjected  to  the  phenomenon. 

Figures  10  through  IS  show  grid  and  axial  velocity  plots  for  the  same  rod 
with  a  smaller  velocity  gradient.  In  this  case  the  gradient  is  1x10*  sec 
i.e.,  the  tail  velocity  is  zero  and  the  tip  velocity  is  1x10*  cm/sec.  Necks 
and  stabilized  regions  are  seen  to  form  in  the  same  manner  as  previously  but 
with  different  sizes  and  velocities.  The  last  three  figures — 130,  140,  and 
150  microseconds — show  in  detail  the  formation,  stabilization,  and  growth  of 
the  second  set  of  elastic  regions. 

Changing  the  length  of  the  rod  has  no  effect  other  than  to  provide  more 
rod  within  which  necks  can  form.  Figure  16  is  a  velocity  versus  axial  posi¬ 
tion  plot  at  150  microseconds  for  a  rod  with  the  same  1x10*  sec~^  velocity 
gradient  but  with  length  increased  to  20  cm.  The  stabilized  velocities  and 
positions  are  identical  to  those  in  Figure  15.  Whereas  the  10-cm  rod  will 
form  perhaps  one  more  set  of  necks  and  a  neck  at  the  center  as  it  continues  to 
stretch,  the  20-cm  rod  will  form  many  more.  Fignre  17  shows  the  neck  system 
for  the  20-cm  rod  at  250  microseconds.  Bight  necks  are  clearly  visible,  and 
more  necks  will  form  as  the  rod  continues  to  stretch. 

The  effect  of  varying  yield  strength  from  2  to  5  kilobars  in  a  staballoy 
rod  is  seen  in  Figures  18  and  19.  These  are  grid  plots  (amplified  in  the 
radial  direction)  at  various  times  for  rods  initially  10-em  long,  1  cm  in 
diameter  with  a  velocity  gradient  of  1x10*  vm/sec.  All  are  assumed  to  have 
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Figure  10.  Grid  and  Velocity  Plot  at  20  Microseconds  -  5-Kb 
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Copper  with  1x10  sec  Gradient 
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Grid  and  Velocity  Plot  at  50  Microseconds 
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Copper  with  1x10  sec  Gradient 
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Figure  12.  Grid  and  Velocity  Plot  at  90  Microseconds 

4  -1 

Copper  with  1x10  sec  Gradient 
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figure  18.  Grid  Plots  for  lO-cm-Long  Uranium  Rods  at 
200  Microseconds 
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I'l.iNSOn's  ratio  of  0.3. 

Fignxe  IS  presents  the  rods  at  200  nicroseconds  for  three  cases.  The  top 
plot  is  for  an  elastic/perfectly  plastic  yield  strength  of  2  Kb,  the  middle 
rod  raises  this  to  5  Kb.  and  the  bottom  plot  is  for  a  linearly  work  hardening 
material.  The  work  hardening  material  initial  yield  strength  is  2  Kb.  and  it 
rises  to  a  saturation  level  of  S  Kb  at  a  strain  of  SO  percent.  The  work 
hardening  rod  and  the  S~Kb  rod  appear  very  much  alike  at  this  time.  This 
similarity  is  expected  since  the  work  hardening  material  is  limited  to  5  Kb  in 
strength.  The  2-Kb  rod  is  significantly  different.  The  elastic  segments  are 
not  as  large  as  those  in  the  other  rods,  and  the  number  of  necks  is  reduced  at 
this  time.  Figure  19  presents  later  time  plots  for  these  same  cases.  Again, 
the  5  Kb  and  the  work  hardening  rods  are  very  similar  at  times  close  to  300 
microseconds.  The  time  in  seconds  is  indicated  on  each  plot.  The  2-Kb  rod  is 
pictured  at  400  microseconds.  It  still  has  a  long  central  section  stretching 
plastically,  whereat  the  other  rods  have  virtually  completed  necking. 

Figure  20  is  the  same  staballoy  rod  at  400  microseconds  bnt  with  no  limit 
on  the  work  hardening  yield  strength.  A  tangent  modulus  of  4  Kb  it  used,  bnt 
there  is  assumed  to  be  no  saturation  stress  level  in  the  material.  Neck 
formation  is  considerably  retarded  over  the  2  Kb  case.  Necks  have  begun  to 
form  at  400  microseconds  but  are  of  very  limited  radial  extent. 

A  xero  yield  strength  copper  rod  stretching  under  a  1x10*  sec  ^  gradient 
is  shown  in  Figures  21  and  22.  There  is  only  slight  velocity  retardation  at 
the  ends,  and  there  is  no  necking  whatsoever.  The  grid  plots  in  Figure  21 
show  the  rod  at  0,  150  and  230  microseoondt.  The  velocity  versus  axial 
position  plot  in  Figure  22  is  st  230  microseconds. 
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Figure  21.  Grid  Plots  for  Copper  Rod  with  Zero  Yield  Strength 


The  calcoletlons  presented  in  Figures  23  end  24  illustrate  the  alaiost 
nonexistent  effect  of  material  sonnd  speed  on  rod  deformation.  Figure  23 
shows  one-half  of  a  copper  rod  with  a  yield  strength  of  5  Kb  and  the  normal 
sound  speed  of  4x10*  cm/sec  at  200  microseconds.  The  velocity  gradient  is 
1x10*  om/sec.  and  the  rod  was  initially  10  cm  in  length.  At  200  microseconds 
just  one-half  of  the  rod  is  almost  10  cm  in  length.  Two  necks  exist,  and  the 
stabilised  elastic  segments  are  travelling  at  3.4x10*  cm/sec,  4.1x10*  cm/sec, 
and  S.OxlO*  cm/ sec.  Elastic  velocity  for  the  three  segments  seen  in  the  half¬ 
rod  picture  are  3.4x10*  cm/sec.  4x10*  cm/ sec,  and  4.85x10*  cm/sec.  These 
velocities  are  within  3  percent  of  those  seen  in  Figure  21.  Clearly  then, 
sonnd  speed  has  little  effect  on  the  rod's  deformation. 

Varying  density  does  have  an  effect  as  one  can  see  by  comparing  the 
staballoy  rod  in  Figure  18  with  the  5-Kb  copper  rod  in  either  Figure  23  or  24. 
It  is  more  dramatically  illustrated  in  Figures  25  through  28.  These  figures 
are  one-half  rod  grid  and  velocity  plots  for  an  aluminum  rod  stretching  under 
a  1x10*  seo~^  gradient.  The  plots  show  the  lower  velocity  end  of  the  rod  at 
40  ,  80.  140,  and  200  microseconds.  The  aluminum  has  a  density  of  2.7  gm/cc,  a 
bulk  sonnd  speed  of  5.4x10*  cm/ sec,  a  shock  velocity/particle  velocity  slope 
of  1.34,  a  Gruneisen  ratio  of  2.1,  a  Poisson's  ratio  of  0.33,  and  a  yield 
strength  of  3  Kb.  At  200  microseconds  a  5-Kb  copper  rod  with  identical 
initial  conditions  has  formed  6  elastic  segments  with  the  first  segment  from 
each  end  being  approximately  5  cm  in  length.  The  3-Kb  aluminum,  at  200 
microseconds,  has  formed  only  2  elastic  segments,  and  each  segment  is  almost  8 
cm  in  length.  The  grid  and  velocity  plot  at  140  microseconds  (Figure  27) 
shows  that  a  second  set  of  elastic  segments  briefly  formed  but  were  not 
sustained. 

In  the  next  section  a  simplified  model  will  be  presented  which  explains 
all  of  the  phenomena  seen  in  these  rod  calculations  with  density  and  yield 
strength  as  the  only  materiel  properties  affecting  deformation  and  necking. 
The  calculations  and  the  model  both  predict  that  oeoling  advances  from  the 


Figure  23.  Grid  and  Velocity  Plot  at  200  Microseconds 
5-Kb  Copper  with  1x10^  sec  ^  Gradient  and 
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Figure  24.  Grid  and  Velocity  Plot  at  200  Microseconds  -  5-Kb 
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Figure  25.  Grid  and  Velocity  Plot  for  3-Kb  Aluminum  at 
40  Microseconds 
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Figure  28.  Grid  and  Velocity  Plot  for  3-KL  Aluminum  at 
200  Microseconds 
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ends  of  the  rod  slowly  toward  the  center.  Observation  of  jets  froa  henis- 
pherical  liners'  Indicates  that  this  aodel  does  duplicate  phenoaenon  seen  in 
initially  constant  diaaeter  rods  with  linear  velocity  gradients.  Hie  calcnla- 
tions  and  aodel  indicate  that  necking  and  fracture  which  initiates  at  any 
other  locations  first  is  caused  by  sharp  changes  in  shape  (aass)  or  in  the 
velocity  gradient  Itself. 
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SECTION  III 
A  SinUFIBD  ■)!«. 


la  thia  section,  a  aodel  will  be  developed  which  e^lains  the  phenoaena 
seen  in  the  rod- stretching  cslcnletions.  Conpsrisoas  with  the  cslcnlations 
indicates  that  the  model  Is  eztrenely  accurate  in  predicting  nany  of  the 
details  of  defornation  and  necking.  It  is  currently  Halted  to  Initially 
constant  diameter  rods  with  linear  velocity  gradients.  The  fnadaaental 
differential  equations  include  all  cases  of  interest,  and  it  is  a  relatively 
simple  task  to  expand  the  aodel  beyond  these  limitations.  The  model  e^lains 
necking  and  breakup  phenomena  observed  in  shaped-charge  jets  with  recourse  to 
none  but  the  most  basic  physics  equations. 

Figure  29  defines  nomenclature  to  be  used  in  aodel  development.  The 
figure  shows  a  velocity  versus  axial  position  plot  at  a  tine  when  nany  necks 
have  been  formed.  Elastic  segments  exist  between  axial  positions  X.  >  and 

IpU 

X.  .  where  i  is  the  number  of  the  segment  beginning  with  unity  at  the  lower 

id* 

velocity  end.  These  values  are  functions  of  time.  Initial  values  will  be 

Indicated  as  X.  .(0).  Velocity  will  always  be  assumed  to  be  initislly  zero  st 
J 

the  lower  velocity  end  of  the  rod,  and  initial  rod  position  will  be  assumed  to 
be  zero  at  this  location.  These  assumptions  cause  no  loss  of  generality  but 
merely  require  that  constants  be  sdded  to  velocities  and  positions  computed 
using  the  model. 

Ay 

The  velocity  gradient,  is  designated  as  X.  Initial  velocity  gradient 

is  K(0).  T  is  the  flow  stress  of  the  rod  and  can  be  a  function  of  strain 
and/or  internal  energy.  The  density  of  the  rod  is  p  and  is  assumed  to  remain 
constant.  The  cross-sectional  area  is  designated  as  A  and  is  A(0)  initially. 

Consider  the  formation  of  the  first  elastic  segments  at  each  end.  The 
velocity  of  the  first  segment  on  the  low  velocity  end  of  the  rod  will  be 
V^(t).  The  velocity  of  the  first  segment  on  t.te  high  velocity  end  of  the  rod 
will  be  V^(t).  Conservation  of  momentum  in  an  Initially  constant  diameter  rod 
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X.  .  ARE  FUNCTIONS  OF  TIME 

1 ,  K 

INITIAL  POSITIONS  VARY  FROM  X.  -(0)  TO  X„  , 

1,0  n , 1 

K(0)  =  —  (t=0) 

0  X 

Y  =*  FLOW  STRESS  OF  THE  ROD 
p  *  DENSITY  OF  THE  ROD 

A(t)  =  ROD  CROSS  SECTIONAL  AREA  AT  TIME  t 
A(t=0)  =  INITIAL  CROSS  SECTIONAL  AREA 


Figure  29.  Model  Nomenclature  and  Definitions 


stretching  under  ■  lineer  gradient  requires  that: 


V^(t)  -  K(0)Xj^  i(0)  -  Vj(t)  (1) 

In  other  words,  the  velocity  gain  at  the  lover  velocity  end  of  the  rod  nnst  be 
eqnal  to  the  velocity  loss  at  the  high  velocity  end. 

There  are  two  aethods  available  for  solving  for  V^(t)  and.  through  Equa¬ 
tion  (1),  for  V  (t).  The  siaplest  is  to  consider  a  balance  of  forces  across 
n 

the  interface  between  the  elastic  end  segments  and  the  plastic,  still- 
stretching  central  section.  Considering  the  lover  velocity,  and  setting 

elastic  segment  deceleration  eqnal  to  the  force  at  the  interface  provides  the 
equation: 

av.(t) 

M  -jJ -  -  YA(t)  (2) 

where  I  is  the  yield  strength  of  the  material  at  the  interface  and  M  is  the 
mass  of  the  constant-velocity,  elastic  segment. 

The  mass,  M,  can  be  moat  easily  written  in  terms  of  the  mass  in  initial, 
or  Lagrangian,  coordinates. 

N  -  pA(0)I^  j(0)  (3) 

Equation  (3)  assumes  that  ^(0)  is  0  and  that  the  rod  has  a  constant  diame¬ 
ter  at  t-0.  The  constant-volume,  plastic  stretching  in  the  central  rod  sec¬ 
tion  requires  that: 

A(t)  -  A(0)/(1+K(0)t)  (4) 

substituting  Equations  (3)  and  (4)  into  Equation  (2)  provides: 
dV  <t) 

pA(0)X^  j(0)  -  -  IA<0)/(l+f(0)t)  (5) 
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Hie  quentity  ^(0)  cen  be  written  as: 

Xl  l(0)  -=  Vj(t)/X(0) 
which  redaces  Eqoation  (S)  to: 


dVi(t) 

Vi(t)  -  =  (J)  K(0)/(1+K(0)t)  («) 


Bqaation  (6)  can  be  integrated  from  t=0  to  t  and  to  and  yields  the 

following  relationship  for  for  the  case  in  which  T  is  constant. 


Vj(t) 


Ln(l+X(0)t) 


1/2 


(7) 


The  velocity,  V  (t),  of  the  upper  end  segaent  will  be  the  tip  velocity  of  the 
tt 

rod  ainus  V^(t).  The  strain  in  the  center  rod  section  will  be  exactly  equal 
to  the  strain  at  the  elastic  interface  until  necking  becoaes  predoainant.  The 
aceuBulated  strain  in  the  center  section  can  be  written  as: 


s  =  -Ln  (A(t) /A(0) ) 


(8) 


With  substitution  of  Equation  (4),  the  strain  can  be  written  as: 

e  =  -Ln  (1/(1+K(0)t))  =  Ln  (1+X(0)t)  (9). 


For  a  work-hardening  aaterial,  then,  the  value  of  Y  at  the  interface  can  be 
written: 


Y  =  Yo  +  T„  Ln  (1+I<0)t)  (10) 

where  Y^  is  the  initial  yield  strength  and  Tj|  ia  the  tangent  aodulus.  Substi¬ 
tution  of  Equation  (10)  into  Equation  (d)  and  integration  yields: 


Vj(t) 


2La  (1+K(0)t) 


y  T 

-^  +  Ln  (l-i-E(O) 
P  2p 


•’ll 


1/2 


(11) 
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If  T  is  also  a  inaction  of  iaternal  energy  density.  Equation  (6)  can  be 
solved  by  including  tbe  appropriate  functional  relationship.  The  internal 
energy  is  siaply  the  strain  energy  developed  as  the  rod  stretches.  It  can  be 
written  as: 

/r(e,I)edT/M  (12) 


where  t  is  the  volume  under  consideration.  N  is  the  aass  of  aaterial  in  the 
voluae,  and  I  Is  the  internal  energy  density.  The  difficulty  of  Incorporating 
energy  dependence  into  Equation  (d)  ia  deteralned  by  the  complexity  of  the 
functional  relationship  between  T  and  I. 

Consider  the  case  of  a  very  siaple  relationship: 

I  -  (1  -  I/IJ  (13) 

O  H 

where  1||  is  melt  energy  density.  The  Internal  energy  density,  I,  is  that 
existing  at  the  interface  which  is  equal  to  that  existing  over  the  nniforaly- 
stretching  center  section  of  the  rod.  For  this  relatively  siaple  case 
Equation  (12)  can  be  written  as: 

I  -  Y  s/p  (14) 

-  I  (1-I/IJ  Ln  (l+K(0)t)/p 
o  a 

Equation  (14)  can  be  solved  for  I  to  provide: 

T^La(l+K(0)t)/p 

^  “  1+I„Y  Ln(l+I(0)t)/p 

a  o 


Equation  (13)  then  becoaes: 

Y^^Ln(l+I(0)t)/p 

^  “  *0  ■  l+I„Y„Ln(l+I(0)t)/p^ 

a  o 


(16) 
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Subttitating  Equation  (16)  into  Equation  (6)  and  iato^xatint  yields: 


Vj(t) 


21, 


Ln(lt^E(0)t) 


—  Ln(l+K(0)t) 

1 - e__ - 


I„(l  +  _  Ln(l+I(0)t)> 

Pl„ 


1/2 


(17) 


Equations  (7).  (11,)  and  (17)  were  coaipared  with  values  of  from  rod 
calculations,  several  of  which  were  discussed  in  the  previous  section.  The 
comparisons  are  contained  in  Table  1.  All  comparisons  are  within  6  percent, 
and  most  are  within  1  percent. 

The  position,  ^2^j(t),  of  the  interface  can  be  predicted  from: 

Xj  j(t)  =  Xj  j(0)(l+I(0)t)  (18) 

where : 


Xl,l(0)  =  Vj(t)/X(0) 

The  accuracy  of  this  position  prediction  is  again  within  a  few  percent  until 
very  late  tines  when  it  diverges  slightly  from  the  calculations.  The  diver¬ 
gence  is  caused  by  the  necking  process.  Equation  (18)  assumes  that  stretching 
continues  at  a  constant  rate.  When  a  second  elastic  segment  is  formed,  the 
comparison  begins  to  diverge  because  the  local  gradient  is  no  longer  exactly 
equal  to  the  gradient  la  the  central,  plastic  region  of  the  rod.  Table  2 
provides  a  comparison  between  some  calculated  I.  .(t)  values  and  those  pre- 
dieted  by  Equation  (18). 
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raut  1.  ■OBBLyCALCDLATlOM  OIWAIISCMS  Rtt 


Material 

Mods  With  Dianeter  - 

'  1  cn  and  Lenath  -  10  cn 

Calcnlation 
V^(kn/ sec) 

I(Kb) 

K(0)(seo"^) 

T(psec) 

Sinple  Model 
V^(kn/see) 

AliiBlniuii 

3 

1x10* 

40 

0.27 

0.27 

80 

0.36 

0.36 

140 

0.44 

0.44 

200 

0.49 

0.49 

Copper 

2 

1*10“ 

160 

0.20 

0.21 

200 

0.22 

0.22 

Copper 

5  with 

IrlO-* 

50 

0.19 

0.20 

linear  tbarmal  str etching 

and  lu^lxlO*  ergs/gn 

m 

100 

0.24 

0.25 

200 

0.28 

0.28 

Copper 

5  with 

lxl0« 

50 

0.20 

0.21 

linear  themal  stretching 

and  I„“2rl0»  ergs/gn 

m 

100 

0.26 

0.27 

150 

0.29 

0.30 

200 

0.31 

0.31 

Copper 

5 

2x10^ 

100 

0.35 

0.35 

Staballojr 

5 

lxl0« 

200 

0.24 

0.24 

322 

0.25 

0.25 

Staballo^ 

lxl0« 

150 

0.21 

0.21 

635 

0.41 

0.40 

Bodn^ with  Dianeter  -  1  cn  and  Lenath  -  20  cm 

Copper 

S 

lxl0« 

240 

0.37 

0.35 
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TilLB  2.  ■NBL/CALOILmON  OOlPAtlSONS  FOB  ^(T) 
go4»  with  Piwter  »  1  <;■  and  Length  -  10  cm 


Material  Y(Eb) 

K(0)(seo"^) 

T(psec) 

Simple  Model 
^1.1^*^ 

Cal eolation 

Xl  i(t) 

AlnBinna  3 

IxlO-* 

40 

3.8 

3.8 

80 

6.5 

6.5 

140 

10.6 

10.1 

200 

14.7 

14.3 

Copper  2 

IxlO-* 

160 

5.0 

5.1 

200 

6.7 

7.8 

Copper  5 

1x10^ 

20 

1.7 

1.7 

140 

7.4 

7.9 

200 

10.5 

10.1 

Copper  5  with 

1x104 

50 

2.9 

3.2 

linear  theraal  softening 

and  I„«lal0» 

ergs/gm 

J1 

100 

4.7 

5.0 

200 

8.3 

8.9 

5  with 

50 

3.0 

3.2 

linear  thermal  softening 

and  I|^2xl0*  ergs/gm 

100 

5,1 

5.2 

ISO 

7.2 

7.8 

200 

9.2 

9.6 

Copper  5 

2x104 

100 

5.3 

5.0 

Staballoy  5 

1x104 

200 

7.2 

8.1 

322 

10.6 

10.9 

Staballoy  Y  =2 

1x104 

ISO 

5.3 

5.3 

T®-6 

M 

635 

30.1 

26.6 

Rods  jrith  Disaiter 

-  1  nnd  XeaxtlL-  20  em 

Copper  5 

1x104 

240 

12.6 

10.2 
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It  is  interesting  to  investigate  the  effect  of  material  properties  on  the 
length  of  the  first  stabilized  segment  at  each  end  of  the  rod.  If  L(t)  is  the 
length  of  this  segment  and  q(0)  is  0,  then: 

L(t)  -  Xi^i(t)  -  Xj  ,j(t) 

=  Xj  j(0)(l+K(0)t)  -  XS*'*** 

y 

-1(0)“*'“”'’  ■  XS"’"' 

Snbstitnting  from  Equation  (7). 

.  f  t  0. 1  a» 

Inspection  of  Equation  (19)  indicates  that  L  decreases  as  T  decreases  and  that 
it  decreases  as  p  increases.  Longer  initial  segments  can  then  be  expected  for 
high  strength,  low  density  materials.  As  a  specific  example,  compare  copper 
with  staballoy.  If  each  has  a  S-Cb-yield  strength,  then  Equation  (19)  pro¬ 
vides  an  L  of  4.S  cm  for  copper  and  3.07  cm  for  staballoy  at  a  strain  of  90 
percent,  assuming  an  initial  gradient  of  1x10*  sec"^.  If  both  have  a  l-ltb 
yield  strength,  then  the  L  for  copper  is  2  cm  and  that  for  staballoy  is  1.4  cm 
at  a  time  when  the  strain  in  the  plastic  section  of  the  rod  is  90  percent. 
The  time  for  a  strain  of  90  percent  can  be  calculated  from  Equation  (9),  given 
the  gradient.  For  E(0)*lxl0*  sec  ^  this  time  is  146  microseconds.  Comparing 
copper  and  aluminum  shovs  that  the  aluminum  length  is  always  greatest.  For 
example,  at  200  microseconds,  a  3-Eb  copper  has  an  L  of  4.9  cm  whereas  a  3-Kb 
aluminum  would  have  an  L  of  7.14  cm,  given  an  initial  gradient  of  1x10*  sec 
If  the  aluminum  also  had  a  5-Kb-yield  strength,  its  stabilized  segment  length 
would  be  (5/3)^^^  times  the  length  at  3  Cb,  or  9.2  cm. 

This  result  is  reasonable  physically.  It  says  simply  that  more  rod  will 
become  elastic  if  the  yield  strength  is  greater  or  if  the  mass  of  the  rod  is 
lower.  Whether  this  is  advantageous  from  a  Jet  penetration  standpoint  is 
another  issue.  In  the  zeal  world  it  appears  that  tradeoffs  occur  betweeu 
density  and  yield  strength.  For  example,  staballoy  is  more  dense  than  copper. 


y 
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but  it  also  has  a  higher  yield  strength.  Copper  is  aore  dense  than  alnainna, 
but  copper  has  a  higher  yield  strength.  Effects  within  a  target  favor  the 
higher  density  aaterial  even  though  necking  and  fracture  aay  be  aore  inhibited 
in  a  lighter  density  aaterial. 

These  equations  can  also  be  used  to  investigate  the  effects  of  velocity 

gradient  on  the  length  of  the  first  elastic  segaents.  The  length  decreases  as 

velocity  gradient  increases.  This  again  is  physically  reasonable  since  the 

relief  wave  velocity  is  opposed  by  the  stretching  velocity  of  the  rod.  For 

example,  the  length  of  the  first  elastic  segaent  is  4.5  ca  for  copper  at  a 

yield  strength  o^  S  Kb  and  a  strain  of  90  percent  and  a  velocity  gradient  of 
-1 


IzlO*  sec  .  The  length  becoaes  2.3  ca  if  the  gradient  increases  to  2x10* 
sec  and  it  drops  to  1.1  ca  if  the  gradient  is  4x10*  sec  For  staballoy 
under  the  saae  conditions,  the  lengths  are  3.1.  1.6  and  0.8  ca. 


It  was  aentioued  earlier  in  this  section  that  there  are  two  aethods  for 
deriving  Equation  (7).  The  first  aethod  was  balancing  force  across  the 
elastic/plastic  interface.  The  second  aethod  involves  balancing  energy  for 
the  entire  rod.  The  aethod  will  be  presented  here  for  the  insight  it  provides 
into  energy  transfer  aechanisas. 

The  initial  energy  in  the  rod  is  given  by: 


E(o) 

Substitution  of: 


J yv»pA(0)dx(0) 


V  -  K(C)x(0) 


provides  the  equation: 


E(0) 


/  "n.: 


(0) 


pA(0)  K(0)*  1(0)*  dl(0) 


(20) 


Integration  yields: 

E(0)  =  J  pA(0)  I»{0)  X*  i(0) 


(21) 
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whore  it  has  been  assnmed  that  ^(0)  =  0,  i.o.,  that  the  initial  Lagrangian 
coordinate  of  the  rod  is  0. 

The  energy  in  Equation  (21)  must  be  balanced,  at  any  tisie,  by  the  kinetic 
and  internal  energy  in  the  rod.  The  kinetic  energy  in  the  still-stretching 
plastic  center  of  the  rod  can  be  written  as: 


KEj  =  pA(0)  K(0)»  -  X}  j(0))/6. 


The  internal  energy  of  the  section  is: 


I  f 
's  ■  A,  . 


Y  Ln  (1+K(0)t)  A(0)  dX(0) 


I,  «  YA(0)  (X,  ^(0)  -  X,  ,{0)  Ln  (1+K(0)t) 
2  n,o  1,1 


The  kinetic  energy  of  the  stabilized  elastic  segment  at  the  lower  velocity  end 
of  the  rod  is: 


KEj  «  j  pA(0)  X,  ,(0)  V? 


1,1""  ’1 


assuming  Xj^  ^(0)  «  0.  The  internal  energy  of  this  segment  can  be  approximated 
by  assuming  that  the  axial  stress  and  strain  vary  linearly  from  zero  values  at 
X(0)»0  to  the  plastic  values  at  X(0)bXj^(0)>Vj^/K(0).  The  stress  at  the 
elastic/plastic  interface  is  the  yield  strength,  T.  The  strain  is 
Ln(l-t-K(0)t).  Therefore: 


h-fo 


o(X)  8(X)  A(0)  dX(0) 


f  ^i,i«»  y: 

•'o  X, 


YX(0) 

X,  ,(0) 


Ln(l^X(0)t)X(0) 

^1,1<0> 


A(0)dX(0) 


YA(0)  Lnd+KCOt)  Xj  j(0)/3. 
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For  the  elaetic  eegaent  at  the  upper  velocity  end  of  the  rod,  the  kinetic 
energy  is: 

IE3  .  |pA(0)  (X^  i(0)  -  l„,o(0))  V*  (26) 

and  the  Internal  energy,  I3  !•: 

I3  -  YA(0)  Ln  (1+K(0)t)  ~ 

setting: 

lEj  +  Ij  +  KEj  +  I2  +  KE3  +  I3  =  E(0) 
and  Making  the  substitutions: 


V 

n 


.(0)  1(0)  -  V, 

Ug  1  1 


and 


n,  o 


*n,l<®^  ~  *1,1 


(0) 


The  equation  can  be  solved  for  as  given  by  Equation  (7).  Energy  is 
balanced  completely  in  Equation  (28)  by  the  appearance  of  the  two  end  elascic 
segments.  It  is  transferred  from  the  upper  velocity  end  of  the  rod  to  the 
lover  velocity  end.  Table  3  provides  internal  and  kinetic  energies  for  the 
three  sections  of  the  rod  at  different  times  for  a  S-Cb  copper  rod,  10  cm  in 
length,  1  csi  in  diameter  with  a  velocity  gradient  of  1x10*  see~^.  The  table 
shows  that  energy  is  transferred  from  the  central  plastically-stretching 
section  to  both  ends  of  the  rod,  providing  the  additional  kinetic  and  internal 
energies  required  to  stabilise  these  end  sections.  The  rate  of  total  energy 
transfer  is  changing  with  time.  The  rate  of  energy  transfer  virtually 
stabilixes  quickly  in  terms  of  energy  density  (energy  per  unit  mass). 
Initially,  the  transfer  rate  is  3.36x10^*  ergs/gm/see.  By  30  microseconds  it 
has  been  reduced,  somewhat,  to  4.8x10*^*,  by  100  microseconds  it  is  reduced  to 
4.66x10**,  and  it  changes  very  little  thereafter. 
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TABU  3.  BOBT  TBAMSnB 


Copper  Rod  vltli  Y  =  S  Kb,  L  =  10  cn.  D  =  1  cm,  K(0)  1x10^  sec 

Total  iaitial  energy  is  1.165x10^^  ergs 


Tine 

(psec) 

(ergs) 

“l 

(ergs) 

(ergs) 

“2 

(ergs) 

^3 

(ergs) 

“3 

(ergs) 

0 

0 

0 

0 

11.65x10** 

0 

0 

SO 

0.1 13x10 « 

0.3 40x10 

0. 913x10 

5.56x10** 

0.113x10** 

4.62x10** 

100 

0.253xl0‘« 

0.7«0xl0‘* 

1.20xl0*« 

4.11x10** 

0.253x10** 

5.07x10** 

150 

0.385xl0« 

1.18.10« 

1 .29xl0*« 

3.26x10** 

0.385.10** 

5.17x10** 

200 

0.505xl0‘» 

1.52xl0‘* 

1.28xl0»* 

2.67x10** 

0.505x10** 

5.17x10** 

250 

0.615xl0‘« 

1. 85x10 « 

1.23x10** 

2.23x10** 

0.615x10** 

5.12x10** 

Tine 

(pseo) 

(ergs) 

I2>KE, 

(ergs) 

(ergs) 

0 

0 

11.65x10** 

0 

50 

0.453x10** 

6.473x10** 

4.733x10** 

100 

1.013x10** 

5.310x10** 

5.323x10** 

150 

1.545x10** 

4.550x10** 

5.555x10** 

200 

2.025x10** 

3.950x10** 

5.675x10** 

250 

2.465x10** 

3.460x10** 

5.735x10** 
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It  it  a  fairly  simple  matter  to  demonstrate  that  energy  remains  balanced 
whether  elastic  segments  exist  beyond  the  end  segments.  The  first  elastic 
segments  on  each  end  totally  aatisfy  all  momentum  and  energy  reqnirementa. 
Were  it  not  for  the  effects  of  neck  corvatnre  on  the  axial  stress  in  the  xoi, 
no  more  necks  or  elastic  segments  wonld  be  formed.  However,  neck  cnrvatnre 
jnst  beyond  the  elastic/plastic  interface  gives  rise  to  a  tensile  hydrostatic 
pressure  which  increases  the  tensile  axial  stress  in  the  neck.  Since  a  stress 
gradient  such  as  this  precipitates  material  acceleration  or  deceleration,  a 
new  wave  is  formed  at  the  neck.  At  the  low  velocity  end  of  the  rod  this 
stress  gradient  is  positive  toward  the  center  of  the  rod,  giving  rise  to  an 
acceleration  in  that  direction.  Conversely,  a  small  deceleration  wave  travels 
into  the  first  elastic  segment.  This  deceleration  has  little  effect  on  the 
relatively  large  first  clastic  segment.  As  it  advances  from  the  neck  toward 
the  center  of  the  rod,  mass  is  accelerated  until  all  of  the  momentnm  in  the 
wave  is  converted.  At  this  point  another  elastic/plastic  interface  is  formed 
and  the  process  is  repeated.  Wave  activity  at  the  upper  velocity  end  of  the 
rod  is  identical  except  that  accelerations  become  decelerations  and  vice 
versa. 


The  necks  at  the  first  elastie/plaatio  interface  begin  to  become  impor¬ 
tant  and  inflnenoe  flow  when  the  acceleration  of  the  interface  drops  below 
that  indnced  by  the  stretching  process.  The  rate  of  change  of  velocity  due  to 
stretching  is  approximated  by: 


=  - 

dt  dx  dt  *  dX  1-M(0)t 

The  neck  region  should  begin  to  develop  curvature  when: 
dV,  V-KO) 


(29) 


dt  1+K(0}t 


(30) 


Differentiating  Sqnation  (7).  enbstitnting  the  derivative  into  Equation  (30) 
and  sinplifying  yields  the  equation: 

[Ln(l+K(0)t)  I  =  2  (31) 

which  indicates  that,  for  elastic,  perfectly  plastic  Materials,  the  time  at 
which  significant  curvature  occurs  is  a  function  of  gradient  only.  The  solu¬ 
tion  to  Equation  (31)  is: 

t  =  (e®*^  -  1)/I(0)  *  0.648/E(0)  (32) 

For  K(0)  =  1x10*  sec  Equation  (32)  predicts  that  neck  curvature  will  begin 
at  t**64.8  microseconds.  Examination  of  rod  calculations  indicates  that  this 
equation  provides  a  reasonably  accurate  estimate  of  the  time  at  which  necking 
becomes  Important,  i.e.,  the  time  at  which  the  axial  stress  begins  to  exceed 
the  flow  stress  of  the  material.  The  equation  is  not  valid  for  work  hardening 
or  thermally  softening  materials.  Equations  for  these  types  of  materials  can 
be  derived  by  combining  equation  Equation  (30)  with  the  appropriate  equation 
for  Vj. 

As  curvature  in  the  first  neck  becomes  important,  the  second  neck  in  from 
each  end  begins  to  form.  The  neck  forms  on  the  central  rod  section  side  of  a 
second  elastic/plastic  front.  The  velocity  of  this  second  elastic/plastic 
front  can  be  found  from  equating  the  acceleration  in  the  first  neck  region  to 
the  difference  in  forces  between  the  first  and  second  elastic/plastic  inter¬ 
faces. 


The  equation  to  solve  is: 


XA(t) 


(33) 


SO 


Snbstitvtiag: 


M  -  p  A(0)  (*2,0^®^  " 

“  P  A<®> 


into  Eqiutloa  (33)  aad  regxoapiag  terms  yields: 


dVi  dVj 

-  Vj)  »  dT  ^  dr> 


2IA(t>I(0)  2Y  K(0) 
pA(0)  “  p  1+I(0)t 


(34) 


Eqaatioa  (29)  cea  be  used  to  approziaste  the  time  derivatives  ia  Bqaatioa 
(34),  redaciag  the  eqnatioa  to  aa  arithmetiic  relatioaship. 

(V2  -  Vj)  (3V^  +  Vj)  =  2I/p  (35) 

The  solatloa  of  Eqaatioa  (35)  is: 

Vj  -  -V^  +  (4V*  +  2Y/p)^''^  (36) 

To  aaiataia  aomeataa  balaaee,  the  correspoadiag  velocity  at  the  upper  ead 
(high  velocity  ead)  of  the  rod  mast  be: 

Vl  “  '<®>»a.l<®)  -  ^2 

After  the  secoad  aeck  begias  to  fora,  a  third  elastic/plastic  iaterface  will 
begia  foraatioa.  The  eqaatioa  for  the  velocity  fox  this  third  iaterface  will 
be : 

V3  -  -Vj  +  (4V*  +  2Y/p)^^*  (38) 

la  geaeral,  thea,  iaterface  velocities  will  be  givea  by  the  relatioaship: 

Vj  -  -Vj_^  +  (4V*_i  +  2Y/p)^^*  (39) 

for  iaterfaoes  propagated  froa  the  rod's  lover  velocity  ead.  These  valves  caa 
be  subtracted  froa  the  iaitial  peak  velocity,  K(0)X  ,(0),  to  obtaia  the 

B,1 

velocity  for  the  correspoadiag  upper-ead  iaterface. 
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Table  4  providea  soae  comparlaoas  between  calculations  and  the  simple 
model  for  these  welocities  at  the  lower  end  of  the  rod.  The  maximum  differ¬ 
ence  is  7  percent,  and  moat  comparisons  are  much  cloaer.  Aa  indicated  in  the 
table,  in  some  cases  segments  have  blended  into  other  segments  as  time  in¬ 
creases. 

Given  a  rod  velocity  gradient  and  velocitiea  for  all  elastic/plastic 
interfaces,  one  can  determine  the  poaition  of  these  interfaces  at  a  given  time 
and  the  mass  between  interfaces.  The  hydrocode  calculations  indicate  that 
these  interface  positions  for  second  and  anbsequent  rod  segments  remain  in  the 
center  of  the  segment  as  each  elastic  region  grows. 

Given  a  set  of  velocity-position  elastic/plastic  interfaces,  one  can 
estimate  a  yield  atrength  and  snbaeqnent  segment  growth  as  well  as  details  of 
the  necking  between  segments. 
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sicriiON  IV 

OCMIGiL  BODS 


The  equations  developed  In  the  previous  section  apply  to  the  sitnation  of 
an  initially  constant  diameter  rod  stretching  under  a  linear  velocity 
gradient.  In  this  section,  the  solution  for  rods  with  a  variable  diameter 
Till  be  explored.  The  specific  case  to  be  addressed  is  that  of  a  rod  which 
initially  is  the  frustum  of  a  cone. 

Figure  30  shows  the  lOOOT  hydrocode  initial  shape  and  zoning  for  a  10-cm 
long  rod  with  an  aft  diameter  of  2  cm  and  a  forward  diameter  of  1  cm.  The 
gradient  is  a  linear  1x10'*  sec  i.e.,  the  forward  tip  of  the  rod  is  travel¬ 
ling  at  1x10*  cm/sec  and  the  rear  of  the  rod  is  initially  at  rest.  Because  of 
the  shape  of  the  rod,  a  velocity  versus  mass  curve  would  not  be  a  linear 
curve.  Instead,  more  mass  would  travel  at  the  lower  velocities.  The  rod  is 
copper  with  a  2-Eb  yield  strength. 

Figures  31  through  38  present  grid  plots  and  axial  velocity  versus  length 
plots  for  the  rod  at  40,  80,  120,  160,  200,  240,  280,  and  300  microseconds. 
As  in  all  previous  cases,  necking  begins  at  the  ends  of  the  rods  and  proceeds 
toward  the  center.  In  this  particular  ease  there  is  not  a  great  deal  of 
difference  between  the  segment  velocities  in  the  calculation  and  those  com¬ 
puted  from  the  simple  equations  of  the  last  section.  For  example,  constant 
diameter  equations  would  predict  that  *  0.25x10*  cm/seo  at  300  micro¬ 
seconds.  The  hydrocode  calculation  predicts  0.224x10*  cm/seo.  The  simple 
equations  predict  V2  *  0.29x10*  cm/ sec  given  the  erroneous  of  0.25x10*,  and 
they  predict  V2  ■■  0.27x105  cm/ sec  given  the  correct  The  hydrocode  calcu¬ 

lation  predicts  that  V2  should  be  OJ6xlO*  cm/ sec. 

It  is  apparent  in  the  calculations  that  the  rod  assumes  a  more  constant 
diameter  shape  as  time  proceeds  so  that  it  would  be  expected  that  elastic 
segment  velocities  would  be  more  nearly  predictable  from  the  constant  diameter 
equations  ae  one  proceeds  away  from  the  ends  of  the  rod. 
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Figure  39  is  the  initial  grid  plot  for  a  more  conically  shaped  rod.  In 
this  case  the  base,  or  aft.  diameter  is  4  cm  and  the  tip  diameter  is  1  om. 
Other  rod  features  are  identical  to  the  previous  case.  Figures  40  through  47 
show  grid  and  axial  velocity  plots  at  40,  80.  120.  160,  200,  240,  280,  and  300 
microseconds.  There  is  far  more  difference  in  this  case  between  the  hydrocode 
calculation  and  the  simplified  constant  diameter  equations.  For  example,  it 
is  seen  that  the  first  neck  has  barely  formed  on  the  low  velocity  end  of  the 
rod  where  two  necks  have  been  formed  on  the  upper  velocity  end.  The  predic¬ 
tion  of  from  the  constant  diameter  equations  at  300  microseconds  is  the 
same  as  in  the  previous  case,  i.e..  0.25x10*  cm/ sec.  The  calculation 

predicts  0.219x10*  cm/sec  on  the  low  velocity  end  and  0.646x10*  ca/sec  on  the 
upper  velocity  end  (instead  of  1x10*  -  0.25x10*  <■  0.75x10*  cm/sec).  Clearly, 
a  4-to-l  diameter  difference  over  the  rod's  length  is  sufficient  to  cause 
large  errors  in  the  predictions  for  the  constant  diameter  equations. 

The  simplified  equations  of  the  last  section  can  be  written  for  variable 
diameter  rods,  but  they  do  not  remain  so  simple.  In  fact,  there  Is  no  obvious 
analytic  solution  to  the  equations,  and  they  must  be  solved  numerically.  In 
this  era  of  programmable  hand  calculators,  this  does  not  present  a  serious 
problem. 

The  volume  of  a  frustum  is  given  by: 

T  -  I*  (*,  * 

where  and  A2  are  the  base  and  top  areas  and  h  is  the  length  of  the  frustum 
along  its  axis.  The  constant  volume  plastic  relation  then  is: 

h(0)  (Aj(0)  +  Aj(0)  +^Ajon^(0))  (40) 

-  h(t)  (A^(T)  +  Aj(T)  +yAj(T)Aj(T)) 
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V 


Now 


MO)  *  i(0)  atsniing  (,(0)  -  0 

tad  h(T)  »  Xi,i(T)  - 

Also,  the  srest  are  faactioas  of  the  laitial  and  final  values  of  X,  ^  and  X.  a 

i»o 

Aj(0)  -  nrj(0)» 

where  r^^CO)  is  the  initial  base  radius. 

(rg(0)  -  rj(0))  2 

Aj(0)  -  wr^CO)*  -  n  rj(0)  +  - ^ -  Xj  j(0)  (41) 

where  r£(0)  is  the  initial  radins  of  the  upper  velocity  end  of  the  rod  and 
is  the  initial  rod  length. 

Since  the  ends  of  the  rod  are  relieved  alaost  instaataneonsly, 
r(T)  -  rj<0) 

and  A^(T)  -  Aj(0)  =  nrj(O)^ 

However,  r2(T)  sinst  be  calculated  froa  Equations  (40)  and  (41)  given  Aj^(O), 
A2(0),  and  A2(T).  The  equation  is  a  quadratic  with  the  only  positive  value  of 
r2(T)  given  by: 

.j(T)  -  -J—  .  1/2  rj(0)*  ♦  <  <;^  -  ,J(0» 

where  the  tern  B  is  given  by: 

B  -  i(O)  (A^(0)  +  A2(0)  +Ja^(0U^0)) 
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lliaa,  the  area  at  the  first  elastic/plastic  interface  is  s  coaplex  fnnction  of 
initial  geoaetry  and  the  position  of  the  interface  in  laboratory  and  Lagran- 
gian  frames.  This  is  the  primary  complication  in  solving  for  an  analytic 
fnnction  providing  the  position  Xj  The  equations  analogous  to  Equations 
(7),  (11).  and  (17)  then  become  very  complex.  There  may  be  substitutions 
available  to  provide  analytic  solutions,  but  they  were  not  found  during  this 
study.  It  remains  then  to  solve  the  equations  numerically.  Numerically 
solving  the  equations  is  a  straightforward  exercise  but  one  which  will  not  be 
pursued  in  this  report. 
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Figure  32.  Grid  and  Velocity  Plot  : 

Dl/D2=2  at  80  Microseco! 


X  C  C  M  ) 


ININ*  I  IH»X«  9  t  JWOI«  tOS 


Figure  33.  Grid  and  Velocity  Plot  for  2-Kb  Copper  Conical  Rod 
D1/D2-2  at  120  Microseconds 
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Figure  37.  Grid  and  Velocity  Plot  for  2-*b  Copper  Conical 
Rod,  D1/D2B2  at  280  Microseconds 
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Figure  38.  Grid  and  Velocity  Plot  for  2 -Kb  Copper  Conical 
Rod,  Dl/D2«2  at  300  Microseconds 
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Figure  40.  Grid  and  Velocity  Plot  for  2-Kb  Copper  Conical 
Rod,  Dl/D2=4  at  40  Microseconds 
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Figure  41.  Grid  and  Velocity  Plot  for  2-Kb  Copper  Conical 
Rod,  Dl/D2*4  at  80  Microseconds 
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Figure  43.  Grid  and  Velocity  Plot  for  2- Kb  Copper  Conical 
Rod,  D1/D2b4  at  160  Microseconds 


Figure  44.  Grid  and  Velocity  Plot  for  2-Kb  Copper  Conical 
Rod,  Dl/D2»4  at  200  Microseconds 
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Figure  45.  Grid  and  Velocity  Plot  for  2- Kb  Copper  Conical 
Rod,  D1/D2b4  at  240  Microseconds 


i igure  47.  Grid  and  Velocity  Plot  for  2- Kb  Copper  Conical 
Rod,  Dl/D2=4  at  300  Microseconds 
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It  is  concluded  froa  this  study  thst  nil  phenoaens  seen  in  the  stretching 
and  deforaation  of  jets  under  large  velocity  gradients  can  be  explained  by 
assuaing  that  the  jet  has  soae  strength,  i.e..  it  is  not  coapletely  aelted. 
It  is  further  concluded  that  careful  ezaaination  of  the  positions  and  veloci¬ 
ties  of  jets  will  allow  predictions  of  strength  to  be  aade  for  very  high 
strain,  high  teaperature  states  of  aetals. 

Further  work  is  indicated  in  the  application  of  the  equations  developed 
in  this  report  to  actual  jets  of  various  aetals.  Soae  predictions  have  been 
aade  for  a  classified  aunition  and  are  included  in  a  separate  letter  report. 
Further  theoretical  work  on  the  equations  theaselves  should  provide  an  under¬ 
standing  of  the  growth  of  the  elastic  portion  of  jet  segaents  beyond  the  end 
segaents.  It  is  believed  that  a  great  deal  of  aaterlal  property  inforaation 
is  contained  within  the  equations  which  predict  this  growth.  There  was  insuf¬ 
ficient  tiae  under  this  project  to  solve  for  those  relationships.  It  would 
also  be  of  interest  to  apply  these  equations  to  stretching  self-forging  frag- 
aents  with  gradients  far  below  those  seen  in  jets. 

This  report  basically  initiates  a  technology  area  which  appears  quite 
fruitful  for  pursuit  in  terns  of  understanding  jet  behavior  and  the  behavior 
of  aaterials  in  very  high  strain  and  teaperature  environaents. 
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